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Abstract 

A procedure is suggested for testing the resolution and comparing the relative accuracy of nu- 
merical schemes for integration of the incompressible Euler equations. 

There is great interest in the behaviour of solutions of the three-dimensional incompressible Euler 
equations [TJ. An important open question is whether the vorticity field can develop a singularity in 
finite time; see, e.g., the review and literature cited in |2J. A central challenge to the computational 
investigation of this question is the inevitable loss of accuracy of any numerical scheme for singular, or 
near singular, solutions. We do not attempt to answer this question ; rather we propose a procedure 
for testing the accuracy and comparing the precision of numerical schemes used to integrate the 
incompressible Euler equations. 

The Euler equations for the velocity field u (divu = 0) are 

— = -Vp, — = d t + u-V, (1) 

Dt Dt ' 

and for the vorticity u> = curl u, 

Dto 

dt u> = curl (u x u>) or = u> • V« . (2) 

The new feature of the proposed test is to introduce a passive tracer concentration 6(x, t) whose 
initial data are under the investigator's control, satisfying 

^=0. (3) 

Dt K ' 

The introduction of 6 allows us to define a variable q = u> ■ V# |3J. An application of Ertel's theorem 
produces the well-known result 

£ = <*> 

Now define 

B=VqxV6 (5) 

which satisfies divB = and comes endowed with initial conditions inherited from those for u and 
9. The fact is (for the proof see the Appendix) that 13 satisfies |4j 

DB 

d t B = curl (it x B) or — = B Vu. (6) 

The vector field B contains information on u, Vw, V# and V 2 # embedded in the combination 
in ([5]). It evolves in the same way as u> and so is subjected to similar (potentially tortuous) stretching 
and folding processes. It can, however, be evaluated at any particular time t in several distinct ways: 
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it can be evaluated from the result of the evolution in © at time t, or it can be computed from 
its definition using u and 9 evolved up to and evaluated at time t. The degree to which these 
distinct evaluations agree or disagree provides a quantitative gauge of the accuracy of the numerical 
computation. It is not clear that there is a natural scale for the inevitable discrepancies produced 
in any particular computation. However, this procedure produces a precise diagnostic quantity that, 
given identical initial data, can be directly compared side-by-side for different numerical computations 
to evaluate their relative accuracy. Here is the suggested test : 

1. Choose initial data for u and 9, thereby fixing initial data for q and B. 

2. Evolve u and simultaneously solve D9/Dt = 0, Dq/Dt = and DB/Dt = B • Vm. 

3. Test the resolution at any time t > by constructing q\(- , t) = u(- , t) ■ V0(- , t) and then : 

(a) compare the solution for £?(•, t) obtained from solving DB/Dt = BVu with B\(- , t) = 
Vq x {- , t) x V0(- , t) 

(b) and, furthermore, compare these with Bi(- , t) = Vg(- , t) x V6?(- , t) where q(- , t) is 
the evolved solution of Dq/Dt = 0. 

4. For fixed initial data for u this procedure may be implemented for a variety of "markers" 
n (- , t) evolving from distinct initial data n {- , 0) to diagnose the numerical accuracy in 
different regions of the flow. 

Because B contains Vo^, comparing the different computations of B, B\ and B2 tests the accuracy 
of the computation of some of the small scale structures in the flow. 

Appendix 

Using conventional vector identities, the proof of ([6j) is (see |4j) : 

Bt = (Vg) t x (V0) + (Vq) x (V0) t 

= -V(u- Vq) x (V0)-(Vg) x [V(u • V0)] 

= - {u ■ V(Vq) + (Vq) ■ Vu+ (Vq) x u} x (V0) 

-(Vg) x {u • V(V0) + (V0) • Vu + (V9) x uj} 
= -u-VB + (Vq)(u-V9)-(V9)(u-Vq) 

+(V0) x (Vg • Vu) - (Vq) x (V9 ■ Vu) 
= curl(uxB) (7) 
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